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ABSTRACT 


One of the more surprising findings after the first year of JWST observations is the large number 
of spatially extended galaxies (ultra-red flattened objects, or UFOs) among the optically-faint galaxy 
population otherwise thought to be compact. Leveraging the depth and survey area of the JADES 
survey, we extend observations of the optically-faint galaxy population to an additional 112 objects, 56 
of which are well-resolved in F444W with effective sizes, Re > 0.25", more than tripling previous UFO 
counts. These galaxies have redshifts around 2 < z < 4, high stellar masses (log(M,/Mo) ~ 10 — 11), 
and star-formation rates around ~ 100 — 1000M5/yr. Surprisingly, UFOs are red across their entire 
extents which spatially resolved analysis of their stellar populations shows is due to large values of 
dust attenuation (typically Ay > 2 mag even at large radii). Morphologically, the majority of our 
UFO sample tends to have low Sérsic indices (n — 1) suggesting these large, massive, optically faint 
galaxies have little contribution from a bulge in F444W. Further, a majority have axis-ratios between 
0.2 < q < 0.4, which Bayesian modeling suggests that their intrinsic shapes are consistent with being 
a mixture of inclined disks and prolate objects with little to no contribution from spheroids. While 
kinematic constraints will be needed to determine the true intrinsic shapes of UFOs, it is clear that an 
unexpected population of large, disky or prolate objects contributes significantly to the population of 


optically faint galaxies. 
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1. INTRODUCTION 


One of the more challenging populations of galaxies to 
characterize over the past few decades has been the op- 
tically dark/faint or *HST-dark" galaxies, where strong 
dust attenuation causes a lack of UV-optical emission 
making them difficult or impossible to detect with fa- 
cilities such as the Hubble Space Telescope (HST) (e.g., 
Barger et al. 1998; Hughes et al. 1998; Coppin et al. 
2006; Elbaz et al. 2011; Williams et al. 2019; Umehata 
et al. 2020; Manning et al. 2022; Barger et al. 2022; 
Gómez-Guijarro et al. 2022; Xiao et al. 2023). Surveys 
in the far-infrared (FIR) through millimeter wavelengths 
with instruments such as Spitzer/IRAC, Herschel, the 
James Clerk Maxwell Telescope (JCMT/SCUBA), and 
the Atacama Large Millimeter/Sub-mm Array (ALMA) 
have gradually revealed this population to be primar- 
ily concentrated around z = 2 (e.g., Barger et al. 1998; 
Hughes et al. 1998; Coppin et al. 2006; Elbaz et al. 2011; 
Wang et al. 2012, 2016; Tadaki et al. 2017; Franco et al. 
2018; Wang et al. 2019). These and other sub-groups of 
optically faint and/or IR/sub-mm detected galaxies fall 
into the broad category of dusty star-forming galaxies 
(DSFGs), whose discovery, classification and properties 
are reviewed in Casey et al. (2014). In addition to pos- 
sessing large quantities of dust, measurements have re- 
vealed moderate to high stellar masses (M. > 10? M) 
and high star-formation rates (SFRs, e.g., Wang et al. 
2019; Elbaz et al. 2011) with a small subset undergo- 
ing intense bursts of star-formation with substantially 
higher SFRs of = 1000/Mg /yr (e.g., Walter et al. 2012; 
Riechers et al. 2013; Marrone et al. 2018; Ma et al. 2019). 
This obscured galaxy population appears to be an im- 
portant, and possibly dominant contributor to the cos- 
mic star-formation rate density out to z ~ 4 (e.g., Barger 
et al. 2012; Madau & Dickinson 2014; Zavala et al. 2021; 
Barrufet et al. 2023). 

Morphological studies of sources with detections at 
both optical/UV and sub-mm wavelengths have revealed 
that the DSFG population is often much more compact 
at longer wavelengths than at shorter wavelengths (e.g., 
Hodge et al. 2016; Oteo et al. 2016; Barro et al. 2016; 
Rujopakarn et al. 2016; Fujimoto et al. 2017; Tadaki 
et al. 2017; Nelson et al. 2019; Gullberg et al. 2019). 
However, this is not always the case, as Sun et al. (2021) 
report two spatially extended DSFGs with sub-mm sizes 
exceeding the IR sizes from Spitzer/IRAC data. These 
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two wavelength regimes, rest-frame optical and FIR, are 
commonly taken to represent the distribution of exist- 
ing stellar mass and of star-formation, respectively, sug- 
gesting that star-formation is primarily centrally con- 
centrated in these objects. This observation can be 
explained if these galaxies are in the process of build- 
ing their central stellar bulges (e.g., Tadaki et al. 2017; 
Nelson et al. 2019; Tadaki et al. 2020), a strong indi- 
cator of future cessation of star-formation (e.g., Franx 
et al. 2008; Bell et al. 2012; Lang et al. 2014; Whitaker 
et al. 2017). Lending further empirical evidence to sup- 
port this physical interpretation, many DSFGs at higher 
redshifts (z ~ 3) have FIR sizes that are consistent 
with the optical sizes of quiescent galaxies at lower red- 
shifts (z ~ 2). The leading hypothesis is thus that the 
DSFG population are the plausible progenitors of mas- 
sive quenched systems at z ~ 2 (e.g., Toft et al. 2014; 
Fujimoto et al. 2017; Tadaki et al. 2017; Suess et al. 
2021), which in turn are likely the progenitors of the 
most massive and quenched galaxies observed in the lo- 
cal Universe. Given the relative importance of optically 
faint galaxies in the context of both the cosmic SFR 
budget and our understanding of massive galaxy evolu- 
tion, it is imperative to have a more complete census 
of their morphologies and stellar populations (e.g., Got- 
tumukkala et al. 2023; Williams et al. 2023a). 

Prior to the launch of the James Webb Space Tele- 
scope (JWST, Gardner et al. 2023), spatially-resolved 
studies of the optically faint population have primarily 
been facilitated by interferometers such as the Submil- 
limeter Array (SMA) or ALMA (e.g., Walter et al. 2016; 
Tadaki et al. 2017; Cowie et al. 2018; Franco et al. 2018; 
Gómez-Guijarro et al. 2022; Chen et al. 2014; Hsu et al. 
2017). Due to their limited field of view and lower sen- 
sitivity, studies have therefore been limited to bright, 
lensed, or previously known galaxies. Thus, while inter- 
ferometers can achieve the necessary spatial resolution 
to characterize galaxies at higher redshifts, their small 
field-of-view makes it challenging to observe or discover 
large samples of galaxies. Conversely, infrared telescopes 
like Spitzer or Herschel can observe many galaxies, but 
their spatial resolution is too poor for morphological 
studies of galaxies. 

All of this has changed with the successful launch 
and commissioning of JWST, whose NIRCam instru- 
ment's large field of view (~ 10 arcmin?) and wave- 
length coverage redder than HST have already enabled 
numerous studies of the optically faint population out 
to 4.4um (Rieke et al. 2023a). Populations of these 


optically faint galaxies have been studied in the Cos- 
mic Evolution Early Release Science (CEERS) survey 
(e.g., Barrufet et al. 2023; Gómez-Guijarro et al. 2023) 
as well as in observations of the SMACS0723 cluster 
(Rodighiero et al. 2023) finding that these galaxies are 
generally at redshifts 2 < z < 7 with high stellar 
masses (M, > 10!?M,;) and high dust attenuations 
(Ay > 2). Breaking down the JWST-observed optically 
faint galaxies observed in CEERS, Pérez-González et al. 
(2023) find that the majority of these sources (7196) are 
DSFGs at 2 « z « 6. However, to better understand 
how these galaxies fit into our overall picture of galaxy 
growth, we need to study their structures to better in- 
fer how these galaxies have assembled. Gómez-Guijarro 
et al. (2023) study the stellar masses and morphologies 
of DSFGs at 3 « z « 7.5 and find a subset of highly- 
attenuated (defined as Ay > 1) galaxies whose main 
difference from their total sample is their ~ 30% smaller 
effective size. 

Given the compact optically faint galaxies reported 
in Gómez-Guijarro et al. (2023) and expectations that 
optically faint populations in general might be com- 
pact, the discovery of optically faint galaxies detected 
in CEERS with apparent sizes greater than 0.25" has 
come as a surprise (Nelson et al. 2023). This visu- 
ally striking galaxy population, dubbed Ultra-red Flat- 
tened Objects (UFOs), are noticeably more elongated 
and redder than the general population in both their 
inner and outer regions. Whereas the expectation may 
have been that these heavily obscured galaxies should 
be more compact than optically bright galaxies lead- 
ing to a larger dust column density (as in e.g. Gómez- 
Guijarro et al. 2023), the UFOs are not more compact 
than other galaxies at the same mass. Given the low 
axis-ratios combined with inferred Sérsic indices close 
to unity, Nelson et al. (2023) suggest that it is most 
likely for the UFO population to be disk-galaxies. In- 
terestingly, in a spatially resolved study of a single pre- 
viously known and spectroscopically confirmed z — 2.38 
optically faint galaxy in the A2744 field, Kokorev et al. 
(2023) are able to characterize the optically faint galaxy 
as being a massive (log(M./Mo) ~ 11.3), highly star- 
forming (200 M5 /yr), dusty, and edge-on spiral galaxy 
with a nearly uniform Ay ~ 4 across its disk. This 
galaxy has the same general characteristics of the UFO 
population described in Nelson et al. (2023) and appears 
consistent with having a disk morphology based on the 
visual identification of spiral structure in F277W imag- 
ing. Although there is some evidence that UFOs are 
highly inclined disks, low Sérsic indices and a projected 
axis ratio distribution skewed toward low values can also 
be indicative of an intrinsically prolate population. Ex- 
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isting samples so far have been too small to perform axis 
ratio modeling in order to statistically determine their 
intrinsic shapes, but here with deeper imaging over a 
wider area we can acquire large enough samples to per- 
form statistically robust axis ratio modeling. 

Looking at the observed color gradients (which can re- 
veal spatial variations in stellar population properties, 
thus reflecting galaxy evolutionary and assembly histo- 
ries) of the UFO population, Nelson et al. (2023) find 
that they are consistent with previously observed neg- 
ative color gradients in galaxies with colors becoming 
bluer as you move out from the center (e.g., Tortora 
et al. 2010; Wuyts et al. 2010; Guo et al. 2011; Szomoru 
et al. 2013; Chan et al. 2016; Mosleh et al. 2017; Suess 
et al. 2019a,b; Miller et al. 2023). However, the UFOs 
differ in that they are still red in their outskirts al- 
beit not as red as in their inner regions. Nelson et al. 
(2023) suggest that this could be driven by large quan- 
tities of dust in these objects. This dust would need to 
be distributed throughout most of a galaxy to account 
for the red colors throughout unless there were other 
spatial gradients in the UFO stellar populations (e.g., 
younger dust-obscured stars in the outskirts and older 
dust-obscured stars centrally). Resolved stellar popula- 
tion modelling is needed to determine why these objects 
are so red out to large radii. 

In order to better understand the intrinsic shapes 
of UFOs and to address what drives their observed 
colors, we present the identification and properties of 
112 optically faint galaxies found in the JWST Ad- 
vanced Deep Extragalactic Survey (JADES) observa- 
tions of both Great Observatories Origins Deep Survey 
(GOODS) fields. Of the 112 optically-faint galaxies, we 
identify 56 as being UFOs and, in this paper, we seek to 
understand the morphologies as well as the stellar popu- 
lation and dust content of these enigmatic objects lever- 
aging the > 5x larger sample along with deeper imaging. 
We describe our data and analysis methods in Section 2 
and present the results of the integrated stellar popula- 
tions in Section 3. In Section 4, we discuss the observed 
color gradients in the context of the stellar population 
and dust content of their inner and outer regions. We 
discuss the inferred projected shapes of the UFOs and 
show the constraints that can be made on their intrinsic 
shapes in Section 5. We discuss the implications of our 
results in Section 6 and we summarize our conclusions 
in Section 7. 

Throughout the paper, we assume the WMAP9 
ACDM cosmology with Qm = 0.2865, Qa = 0.7135 and 
HO = 69.32 kms ! Mpc^! (Hinshaw et al. 2013). All 
magnitudes in this paper are expressed in the AB system 
(Oke 1974). 


2. DATA ANALYSIS 
2.1. Observations 


JADES! is a joint collaboration between the NIR- 
Cam and NIRSpec science teams carrying out observa- 
tions of both GOODS fields with NIRCam and MIRI 
imaging as well as NIRSpec spectroscopy (Eisenstein 
et al. 2023). Observations were conducted in 9 differ- 
ent filters (F090W, F115W, F150W, F200W, F277W, 
F335M, F356W, F410M, and F444W) spanning 0.90um 
to 4.4um and reaching flux limits of ~ 30 AB magni- 
tudes in F090W, F115W, F150W, and F200W in the 
GOODS-S field whose data is mildly deeper than in 
GOODS-N (~ 1— 1.5 mags deeper than CEERS). Addi- 
tional NIRCam medium-band data from the JWST Ex- 
tragalactic Medium-band Survey (JEMS; program 1963, 
Williams et al. 2023b) is included. The areas covered by 
the GOODS-S and GOODS-N mosaics are ~ 67 and 
~ 58 square arcminutes, respectively (~ 2x wider than 
CEERS). 


2.2. Data Reduction and Photometry 


We follow methods outlined in the first JADES data 
release (Rieke et al. 2023b) for the data reduction and 
photometry and briefly summarize the main steps here. 
The raw data was processed through the JWST Cali- 
bration Pipeline (v1.8.1, Bushouse et al. (2022)) using 
the CRDS pipeline mapping (pmap) context 1009 . The 
first two stages of the pipeline are followed using the 
default parameters to perform detector-level corrections 
as well as flat-fielding and the flux calibration. Before 
combining individual exposures, several custom correc- 
tions were performed to take care of features associated 
with the NIRCam images. Namely, these are a back- 
ground subtraction using the PHOTUTILS Background2D 
class, removal of 1/f noise associated with image read- 
out, and subtracting the “wisp” features using stacked 
wisp templates from JADES and other programs. Fol- 
lowing astrometric alignment using a custom version of 
JWST TweakReg, the images for a given filter / visit 
are combined with Stage 3 of the JWST pipeline and 
these visit-level mosaics are then combined to produce 
the final mosaic (0.031" /pixel). 

Sources are detected following the methods in Rieke 
et al. (2023b) using a signal-to-noise ratio (S/N) image 
generated from inverse-variance weighted stacks of the 
F277W, F335M, F356W, F410M, F444W science and er- 
ror images as the signal and noise images, respectively. 


1 Some of the data presented in this paper can be found in the 
Mikulski Archive for Space Telescopes (MAST), 10.17909/8tdj- 
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Photometry is performed for all of the sources in the cat- 
alog in circular apertures of diameters 0.2", 0.3", 0.5", 
0.6", and an aperture enclosing 8096 of the total energy. 
Errors are estimated by placing 100,000 random aper- 
tures in groups of 1000 across each band's mosaic and 
measuring the rms flux (in electrons) as a function of 
aperture size which is used to assess the contribution of 
the sky background to the total flux uncertainty. 

Our targets are selected using aperture-corrected, 
PSF-matched photometry with a circular aperture of 
diameter 0.3" and as such we infer the integrated stellar 
populations using the same photometry. For the pur- 
pose of determining the colors and the stellar popula- 
tions of the inner and outer regions of each galaxy, we 
make use of the same inner aperture (without the aper- 
ture correction) to describe the inner region and use an 
outer annulus of width 0.15" (determined by subtracting 
the non-aperture corrected photometry within the in- 
ner 0.15" radius from the outer 0.3" radius) to describe 
the outer region. Additional Hubble Space Telescope 
(HST) photometry from the reductions of the Hubble 
Legacy Fields (HLF, Illingworth et al. 2016; Whitaker 
et al. 2019) is utilized in the SED fitting as well. 


2.3. Sample Selection 


For this work, our goal is to leverage the larger area of 
the JADES survey to generate a larger sample of UFOs 
in an attempt to address their nature (e.g., are they 
disks?). To this end, we select the galaxies in JADES 
that are brightest at the reddest wavelengths of our cov- 
erage (4.4um) and faint at HST wavelengths (< 1.6m). 
Examples of HST images for a selection of our targets 
can be seen in Figure 3. From this sample of optically- 
faint galaxies, we take those with semi-major axis effec- 
tive (half-light) sizes greater than 0.25” as our sample 
of UFOs. 

Originally, we adopted the optically-faint galaxy 
(OFG) selection criteria in Nelson et al. (2023) requir- 
ing AB magnitudes of F444W < 24.5, F115W > 27, 
and F150W > 25.5. However, during our analysis we 
noted that some of the selected UFOs showed blue ex- 
cesses either in the form of bluer clumps, excess blue 
emission on the outskirts, or SEDs more indicative of 
bluer stellar populations. These bluer outliers had in- 
ferred Ay values much lower (< 1) than their redder 
companions, and thus to select a more intrinsically red 
parent sample from which to find UFOs, we revised the 
OFG selection criteria of Nelson et al. (2023) to include 
the same F444W brightness cut, but with the addition 
of the following color cuts: F150W - F200W > 0.75 and 
F200W - F444W > 2.0. The selection change results in 
a much more robust sample of thoroughly red objects 
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Figure 1. F200W-F444W color versus the semi-major axis 
effective size illustrating our optically-faint galaxy color se- 
lection and our UFO size cut (dashed black line at 0.25"). 


whose optical faintness is being largely driven by dust 
attenuation (see Section 3). 

With this criteria, we identify 112 total galaxies with 
a range of apparent sizes. We fit for their morpholog- 
ical parameters (described in Section 2.5) allowing us 
to identify UFOs amongst our sample of optically-faint 
galaxies. This selection identifies 56 UFOs whose stellar 
populations, morphologies, and comparisons with our 
broader optically-faint sample we discuss below. In Fig- 
ure 2 we show F200W, F277W, and F444W color images 
of the UFOs and in Figure 3 we show, for a sub-sample of 
UFOs, how HST color images (F606W, F125W, F160W) 
compares with the corresponding JWST image clearly 
illustrating the ability of JWST to reveal galaxies that 
were invisible or faint at HST wavelengths. 

Additionally, from the F444W < 24.5 mag parent 
sample, we identify those that are bluer than our op- 
tically faint selected sample in both colors and analyze 
them as well. This sample is further limited to the same 
redshift and stellar mass range as the whole optically- 
faint sample (described in Section 3). We fit for the 
stellar populations and morphologies of this sample in 
the same way as for the optically-faint galaxies in order 
to more broadly compare the two populations. Unless 
otherwise stated, figures will refer to the blue sample of 
F444W-selected galaxies as the F444W-parent sample 
and the smaller (non-UFO) optically-faint galaxies as 
OFGs. Our three samples are shown in Figure 1 where 
we plot their F200W - F444W colors versus their semi- 
major axis-effective size illustrating both our size cuts 
and one of our color cuts. 


2.4. Stellar Populations 


We utilize the photometric redshift code, EAzY 
(Brammer et al. 2008), to infer the redshifts of our sam- 
ple of galaxies. The templates and photometry used to 
constrain the redshifts follow the methods outlined in 
Hainline et al. (2023). For the rest of the stellar pop- 
ulation parameters, we use the PROSPECTOR SED fit- 
ting code (Johnson et al. 2021) fixing the redshift to 
the value inferred from EAzY. We verified that allow- 
ing the redshift to remain free in our PROSPECTOR fits 
yields redshifts consistent with EAzY as well as yield- 
ing consistent stellar population parameter inferences as 
the fixed-redshift results. To vastly decrease the time 
needed to perform each SED fit, we make use of a neu- 
ral net emulator, PARROT, that has been trained to infer 
the photometry for some underlying stellar population 
model (see Mathews et al. (2023)). This emulator gives 
consistent results with standard PROSPECTOR, but with 
a factor of ~ 10? increase in speed. 

We adopt the PROSPECTOR-{ model described in de- 
tail in Wang et al. (2023) adjusted to keep the redshift 
fixed. In brief, this is the same model as the PRosPEC- 
TOR-o model (Leja et al. 2017), but makes use of a joint 
prior on redshift, stellar mass, and stellar metallicity dis- 
favoring high mass, high-z solutions while still allowing 
a non-negligible probability of obtaining such solutions 
(thus making it possible to discover massive galaxies 
at high redshifts). The stellar populations are modeled 
using the MILES spectral templates (Sánchez-Blázquez 
et al. 2006; Falcón-Barroso et al. 2011) with the MIST 
stellar isochrones (Choi et al. 2016) as implemented in 
FSPS (Conroy et al. 2009, 2010). We adopt the Chabrier 
initial-mass function (IMF, Chabrier (2003)) and model 
the star-formation history (SFH) by fitting for the mass 
formed in seven logarithmically-spaced time bins us- 
ing a continuity prior (Leja et al. 2019a) that weights 
against sharp changes in the SFR between adjacent time 
bins. In practice, this prior is a Student's-t distribu- 
tion for the log ratio of the SFR in adjacent time bins, 
log(SFR, /SFR,..1), with width of ø = 0.3, and degrees 
of freedom, v — 2. 

To accurately characterize the red and likely dusty na- 
ture of these galaxies, we adopt the two-component dust 
model of Charlot & Fall (2000) which characterizes the 
dust in galaxies with a diffuse dust component attenu- 
ating all stars equally and a separate birth-cloud com- 
ponent providing additional dust attenuation for stars 
younger than 10 Myr. Each component can have ex- 
tinctions up to Ay ~ 4.3 and we also allow for a varying 
dust attenuation curve following the prescription in Noll 
et al. (2009). We additionally test two alternative dust 
models in which we: 1) remove the young star compo- 
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Figure 2. Color images (F200W, F277W, F444W) of 35 of our UFOs illustrating the extended red nature of these objects. 
Identifiers in the upper left are the IDs in the JADES photometric catalogues and a 1" bar is shown for reference in the 
lower right. These 1" scale-bars translate to physical sizes of 6.7 kpc to 8.5 kpc depending on the redshift of the object. For 
comparison, we show a "typical" galaxy in the upper left panel with bluer colors and more pronounced color variations than 
the UFOs. 


Figure 3. Comparison between the HST color image (F606W, F125W, F160W) and the JWST color image (F200W, F277W, 
F444W) for a representative sample of our UFOs. As can be seen, these galaxies are very bright and prominent with JWST 
imaging, but almost or completely invisible at HST wavelengths. 
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nent, and 2) fix the slope of the attenuation curve to the 
Calzetti et al. (2000) value. These two variations give 
consistent results with our fiducial dust model. The ef- 
fects of dust emission are included using the dust emis- 
sion templates from Draine & Li (2007) with three free 
parameters that control the shape of the infrared (IR) 
SED (see Leja et al. (2017) for details). Finally, the 
possibility of an AGN powering hot dust emission is in- 
cluded (Leja et al. 2018) utilizing the AGN templates of 
Nenkova et al. (2008a,b) and nebular line emission is in- 
cluded following the CLOUDY (Ferland et al. 1998, 2017) 
FSPS implementation described in Byler et al. (2017). 
We note that certain aspects of this model (i.e., far- 
infrared dust and AGN emission) are not strictly nec- 
essary, but are included as default parameters in the 
emulator we use. It is possible to fix certain parameters 
to narrow the explored parameter space, but doing so 
has negligible impact on our results. 

We fit each object using the photometry as described 
in Section 2.2 and values are reported as the median of 
the posteriors and 1-c uncertainties as the the 84th—50th 
and 50th-16th interquartile ranges. To approximately 
account for systematic errors or issues with early JWST 
photometric calibrations, we enforce an error floor of 5%. 
We show the results of our fits in Figure 4 for a represen- 
tative group of galaxies in our sample. Each fit shows 
the observed SED, the model SED, and a spectrum gen- 
erated in FSPS using the inferred model parameters. 


2.5. Morphologies 


We fit for the semi-major axis effective (half-light) 
sizes and other morphological parameters (Sérsic in- 
dices, n, and axis-ratios, q — b/a) for the entire sam- 
ple of optically-faint galaxies with GALFIT (Peng et al. 
2002, 2010) and with LENSTRONOMY(Birrer & Amara 
2018; Birrer et al. 2021). The setup we use for GAL- 
FIT follows that of Suess et al. (2022) using empiri- 
cal point-spread functions (ePSFs) generated with the 
EPSFBUILDER class in PHOTUTILS (Ji et al. 2023; An- 
derson & King 2000; Anderson 2016). From the mosaics, 
we create 80x80 pixel cutouts centered on each object 
which sufficiently captures all of the light from each tar- 
get. Sources and masks in each cutout are defined using 
PHOTUTILS with any object within 3" of the target cen- 
ter and no more than 2.5 mag fainter than the target 
being fitted. Finally, we subtract the background from 
each image using the SExtractor algorithm in PHOTU- 
TILS and fit for the morphological parameters in both 
the F200W and F444W filters. 

LENSTRONOMY is designed for the modeling of grav- 
itational lenses, but is also equipped to perform image 
modeling in a similar way as GALFIT (i.e., via Sérsic pro- 


file fitting taking into account the PSF). By fitting with 
both LENSTRONOMY and GALFIT, we obtain a rough 
way to validate the morphological fits with the two fits 
tending to give consistent results with few exceptions. 
'The two codes return similar results for the three mor- 
phological parameters (n, q, re), with median differences 
between the two fits of less than 0.1 for each parame- 
ter. Additionally, as discussed in Section 6, we perform 
modeling of the observed axis-ratio distribution to con- 
strain the UFOs intrinsic shapes and by using two sets of 
results we can obtain a rough approximation of the me- 
dian uncertainty on the axis-ratio measurements. Given 
that LENSTRONOMY is able to fit the entire sample of 
UFOs and OFGs whereas GALFIT fails on multiple ob- 
jects, we adopt the LENSTRONOMY results as our fiducial 
morphological parameter values. 


3. INTEGRATED STELLAR POPULATIONS 


In Figure 5, we show how our galaxies are distributed 
in stellar mass and redshift as well as how they are sit- 
uated within the “star-forming main sequence" (SFMS) 
of Leja et al. (2022) for z = 2.5 galaxies. The UFOs 
tend to occupy the same region of both the z — M, and 
SFR — M, plane as their optically faint parent sam- 
ple suggesting that the primary distinguishing factor is 
their structure. T'hese objects tend to be clustered be- 
tween 2 « z « 4, with moderate to high stellar masses 
(log(M./Mo) > 10), and SFRs that tend to lie on or 
above the SFMS of z = 2.5 galaxies from Leja et al. 
(2022) except at the massive end where there is a no- 
ticeable drop to below the SFMS at a mass of ~ 101? M. 
We discuss this drop more in Section 6. Median values 
and percentiles (16th and 84th) for these three param- 
eters (z,log(M./Mo),log(SFR/(Mo/yr)) are 2.42*05, 
10.044963, and 1.581039 for the UFOs and 2.78* 661, 
10.29+9:39, and Las for the smaller optically faint 
galaxies. We show these values along with the values for 
other populations in Table 1. 

Fig. 5 also compares the optically faint populations 
in the context of the bluer F444W-selected parent sam- 
ple and a FIR-selected sample from Ma et al. (2019). 
The FIR-selected galaxies are much more extreme than 
the UFOs with stellar masses predominately higher than 
M, ~ 10!!Mo and SFRs > 100M,/year, and z > 3, al- 
though we note that there is likely a systematic offset 
in their measured SFRs towards higher values compared 
to those we infer here owing to different methodologies 
used (e.g., Leja et al. 2019b). These galaxies are exam- 
ples of extremely red objects with large Ay that were 
selected to have rising Herschel/SPIRE flux densities 
(Soo > $350 > S250). Compared to the bluer F444W- 
selected galaxies, the optically-faint galaxies tend to 
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Figure 4. Observed and model SEDs for a representative sample of our UFOs with the observed JWST photometry as black 
points, observed HST/Spitzer photometry as gray points, model photometry as red points, and the model spectrum in red. 
Photometric points with negative fluxes are indicated with the gray arrows. The SEDs in each column show the lowest (top 
SED) and highest (bottom SED) value galaxy in terms of stellar mass (left), sSFR (middle), and diffuse dust Ay (right). We 
note that some of the SED fits are not ideal, particularly for the HST/Spitzer photometry of the galaxies in the left column. 
However, fits like these are not representative of the majority of our sample, and we also find that the stellar population 
parameter inferences are not strongly changed if only JWST photometry is used, so we keep these objects in our analysis 


have slightly higher redshifts and stellar masses than the 
optically-faint galaxies. The SFMS in Leja et al. (2022) 
is fit to the 1.6um-selected 3D-HST galaxies, so the fact 
that the majority of our sample lies on and around this 
line means that these are not extreme objects but rather 
a normal subset of the massive galaxy population that 
was previously missed. Thus, the population of UFOs 
studied here are not extreme objects in regard to stellar 
mass, star formation rate, or redshift. 

We also show the position of the three FA44W-selected 
samples in the mass-weighted-age (MWA)— Ay plane 
and the stellar-mass Ay plane in Figure 6 with the 
MWA and Ay (for the diffuse dust component) repre- 
senting two likely drivers of optical faintness. Most of 
the optically faint galaxies have large Ay (2.764935 for 
UFOs and 2.61+9:32 for the smaller OFGs) and a young 
MWA (0.9779; for UFOs and 0.901139 for OFGs). 
The right panel additionally shows the inferred Ay in 
the inner and outer regions as discussed in Section 4 
demonstrating the large Ay values in both the inner 


and outer regions of the UFOs. Compared to the bluer 
F444W-selected galaxies, the UFOs have much higher 
dust attenuation, but similar mass-weighted-ages sug- 
gesting that their red colors — and the reason they were 
previously missed is owing to dust attenuation and red- 
shift as opposed to old stellar ages. These results are 
discussed further in Section 6. 

'To further quantify the differences between these pop- 
ulations (UFOs, smaller OFGs, and F444w-parent) in 
these 2D parameter spaces, we perform 2D Kolmogrov 
Smirnov (KS) tests on these distributions, using the 
public code NDTEST?, based on the algorithm described 
in Peacock (1983) and Fasano & Franceschini (1987). 
First, in comparing the F444W-parent sample with ei- 
ther optically-faint population, we find very small p- 
values (< 1074) showing that the optically-faint galaxies 
are a distinct population in terms of stellar mass, SFR, 


? https:/ /github.com/syrte/ndtest 
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dust attenuation, and age. Comparing the UFOs and 
the smaller OFGs, we find a large p-value (0.53) for the 
MWaA-Ay distribution, a moderate p-value (0.14) for 
the log(SFR)-log(M/Mo), and small p-values (< 0.06) 
for the z-log(M/Mo) and Ay-log(M/Mo) distributions. 


This suggests that UFOs are most similar to their 
smaller counterparts in their ages, dust attenuation, and 
SFRs, while stellar mass is where UFOs differ most from 
the smaller OFGs. 


Table 1. Integrated Stellar Population Parameters 


Parameter UFOs OFGs F444W-parent Herschel Ultrared 

z pa 2780 2278 —— 3B 

log(M./Mo) 10.04+9:$3 10.291956 9.811953 [tous 

log(SFR/(Mo/yr))  158'03; 1-481057  — 0941057 2.861517 

by 219135 2.614982 ort ] 

MWA [Gyr] 0.97122 — 0.907120 1.09* 1-55 - 
NoTE— Reported quantity is the median and percentiles (16th and 84th) of the 


distribution for that parameter. 


4. RADIAL PATTERNS IN COLORS AND 
STELLAR POPULATIONS 


Looking at the images of our galaxies, it is apparent 
that these galaxies are consistently red throughout their 
extended light profiles. To make this visual assessment 
more quantitative, we calculate the F200W-F444W col- 
ors for the UFOs in the inner regions (0.3" diameter cir- 
cular aperture) and outer regions (0.6" - 0.3" diameter 
apertures) of the PSF-matched images. We additionally 
calculate the integrated colors for the smaller optically 
faint and bluer F444W-selected galaxies. We show these 
colors in Figure 7 with the UFOs inner and outer colors 
shown with filled and unfilled red circles, respectively. 
As can be seen, the UFOs are strikingly red in both their 
inner and outer regions with both regions colors much 
redder than the bulk of FA44W-selected galaxies. These 
are the same trends that were observed in the sample of 
UFOs first reported in Nelson et al. (2023) where they 
saw that their 12 galaxies had red colors throughout the 
inner and outer regions that were generally much redder 
than the parent population. This is a unique population 
in that the inner and outer regions are both truly red as 
opposed to many galaxies that are observed to have red 
centers and blue outskirts. However, the color gradient 
is still negative in both cases with the UFOs being less 
red in the outskirts than they are in the centers. 

'To further elucidate the persistent redness shown in 
Figure 7, we also infer the UFO stellar populations 
within the inner 0.3" diameter (2.42 kpc at z — 2.5) 
and the outer 0.6" - 0.3" annulus. This approximate 
spatially resolved analysis serves to place constraints on 


what is driving the red colors observed throughout the 
extent of these objects. A representative best-fit spec- 
trum from each region is shown in the lower panel of 
Figure 8. To generate these region-specific spectra, we 
first generate a rest-frame spectrum in FSPS for each 
UFO/region taken from the maximum of the parameter 
posteriors for that objects fit. Then, for each region, 
we stack the individual rest-frame spectra by taking a 
median at each wavelength/flux value to generate a ^me- 
dian" spectrum for each region. It is clear that the two 
spectra are most similar at long wavelengths and least 
similar at shorter/bluer wavelengths with more blue flux 
in the outer regions consistent with the slight color gra- 
dients observed in UFOs that we showed in Figure 7. 

Turning to what drives these slight color gradients, in 
the upper panels of Figure 8, we show the comparison 
between the inner and outer SED fitting results for the 
following parameters: stellar mass (M.), visible atten- 
uation from dust(Ay), mass-weighted age (MWA), and 
specific star-formation rate (sSSFR). All comparisons are 
shown as the ratio or difference of the inner parameter 
inference to the outer parameter inference as a function 
of the stellar mass inferred from the integrated fits (Sec- 
tion 3). The red dashed line in each panel reflects the 
median ratio or difference between the inner and outer 
parameters, the gray dashed line is where the two re- 
gions have the same value, and the red shaded region 
shows the ~ lo spread in the ratio/difference. This in- 
formation is additionally summarized in Table 2, where 
we show the median and percentiles of the inner and 
outer stellar population parameters. 

The largest radial differences are observed between the 
inner and outer Ay, where it is clear that the majority of 
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Figure 5. Left panel: Stellar mass inferred from SED fitting versus the redshift inferred from EaZY with UFOs in red, smaller 
optically-faint galaxies in black, the bluer F444W-selected sample in gray, and ultrared SMGs from Ma et al. (2019) in purple. 
Right panel: how our objects are distributed in SFR and log(M/Mo) with the shaded gray region showing the SFMS at z = 2.5 
from Leja et al. (2022). The UFOs are a much less extreme population than the optically-faint FIR-selected galaxies, with lower 
stellar masses and SFRs. They occupy the high mass end of all FA44W-selected galaxies showing that it is not only extreme 
star-bursts that were missed in previous censuses, but also more moderate galaxies. 


UFOs have higher central Ay than outer Ay although 
both regions have high Ay (Figure 6). Taken at face 
value, this implies that dust concentrations are higher 
in the center as has been found before in studies of other 
galaxy populations (e.g., Nelson et al. 2016; Whitaker 
et al. 2017; Tacchella et al. 2018). On the other hand, 
the gradients in the other parameters, MWA and sSFR, 
are not quite as strong, especially for the MWA which 
has a median ratio of — one, but with significant scatter 
about the median (typical spread of — 0.7). Here, the 
slight color gradients of our UFOs appear to be primarily 
caused by stronger central dust attenuation consistent 
with what was found in Miller et al. (2023) for HST- 
selected star-forming galaxies. 


Table 2. Inner and Outer Stellar Population 


Parameters 
Parameter Inner Outer 
log(M./Mo) 9.86* 025 10117022 
log(sSFR/yr-!) | —8.28*075 | —8.64*0 38 
+1.46 +0.32 
Av 4.26* 120 2.29*0 27 
MWA [Gyr] 0.9470 29 108 Fer 


Table 2 continued 


Table 2 (continued) 


Parameter Inner Outer 


NoTE— Reported quantity is the median and 
percentiles (16th and 84th) of the distribution for 
that parameter. 


5. UFO STRUCTURE 
5.1. Sizes and Sérsic indices 


Inferred morphological parameters are shown in Fig- 
ure 9 and the median/percentiles of these parameters 
are shown in Table 3. In the leftmost panel of Figure 9, 
we show how the UFOs and their optically faint and 
F444W selected parent populations are situated in the 
size-mass plane. For comparison, we show the z = 2 
size-mass relations from Suess et al. (2019a) for star- 
forming galaxies and quiescent galaxies. Our sample of 
UFOs lie on or above the star-forming galaxy size-mass 
relation, showing that these are truly extended objects, 
both physically and apparently. Further, the full sample 
of optically faint galaxies is fairly evenly distributed in 
the size-mass plane relative to the parent sample. This 
is surprising. While we expected optically faint galaxies 
to have compact sizes to drive their high values of dust 
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Figure 6. Left: mass-weighted age (MWA) versus dust attenuation for the UFOs (red), smaller optically faint galaxies (black), 
and the bluer F444W-selected galaxies in gray. The optically faint galaxies have larger dust attenuations than the bluer 
comparison sample, but similar MWAs suggesting that the redness of the optically faint galaxies is largely driven by elevated 
dust attenuation. Right: stellar mass versus dust attenuation for the same samples along with the values inferred for the inner 


and outer regions of the UFOs as described in Section 4. 
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Figure 7. Observed colors versus observed AB magnitude 
in F444W for our optically faint galaxies and the bluer, 
mass and redshift matched, F444W-selected sample. For the 
UFOs, we show the color in the inner 0.3" circular aperture 
(filled red circles) as well as in an outer annulus of width 0.3" 
(open red circles). The optically-faint galaxies are clearly 
much redder than the parent sample with UFOs being con- 
sistently red throughout with red colors in their inner and 
outer regions. 


attenuation, this does not appear to be the case. The 
distribution of the sizes of the optically faint population 
is similar to that of the parent population. 

'The middle two panels of Figure 9 show the distribu- 
tion of Sérsic indices versus stellar mass and effective 
size for the UFOs in the context of all optically faint 
galaxies and the F444W selected parent sample. The 
majority of the UFOs and the optically faint galaxy pop- 
ulation in general have n « 2, suggesting they have 
minimal structural contribution from a bulge. Con- 
trary to our physical expectation that optically faint 
galaxies would be centrally concentrated to drive their 
high dust attenuations, the bulk of these objects have 
Sérsic indices consistent with surface brightness profiles 
that are only slightly steeper than exponential. In- 
terestingly, amongst F444W selected massive galaxies 
(log(M../Mo) > 10) with low Sérsic indices (n < 2), 
and re > 2 kpc, 69% are optically faint UFOs. Thus, 
neither compact sizes nor high Sérsic indices can explain 
the high values of dust attenuation in these galaxies. 
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Figure 8. How the inferred stellar population parameters and predicted spectra of the UFOs compare between their inner 
regions (0.3” circular aperture) and outer regions (0.6” - 0.3” annulus). Top Panels: From left-to-right we show the ratio or 
difference between the inner and outer stellar mass, Av, MWA, and sSFR as a function of the stellar mass inferred from the 
integrated fits discussed in Section 3. The red and gray dashed lines show the median value of the ratio and where the values 
are equal. The red stripe shows the ~ lo spread in the ratio or difference. Bottom Panel: median posterior spectra for the 
inner and outer regions generated from the maximum probability samples of the individual objects showing that the inner and 
outer median UFO spectra are most similar at redder wavelengths. 


Table 3. Morphological Parameter Sum- Table 3 (continued) 
Mary Parameter UFOs OFGs 
Axis-ratio (q) Q.30 019 0.601633 
aided UFOS DEOS Sérsic Index (n) — 112*193 — i,72*198 
log(re/kpc) (dps [ur esed 


Table 3 continued 
Table 3 continued 
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Figure 9. From left to right, we show how the UFOs, smaller optically-faint galaxies (OFGs), and bluer FA44W-parent samples 
are distributed in the planes of size — mass, Sérsic index — mass, Sérsic index — size, and axis ratio — size. The histograms for 
the OFGs include the contribution from UFOs and the FA44W-parent histograms include all galaxies in our samples. The solid 
black lines in the size-mass panel show the size-mass relationships from Suess et al. (2019a) demonstrating that UFOs generally 
lie on or slightly above the star-forming relation while the OFGs are scattered around the quiescent size-mass relation (which 
is largely a result of our size cut). From the other two panels, we can see that the optically-faint galaxies have a much higher 
density of n « 2 objects, particularly for UFOs, relative to the bluer FA44W-parent sample. There are no strong trends between 
the Sérsic index and stellar mass or effective size. The right most panel shows q versus re illustrating that optically-faint galaxies 
have a range of sizes and tend to have lower q values (especially for the UFOs). 
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Figure 10. Axis-ratio distributions for the full F444W- 
parent sample, the optically-faint sample, and the subset of 
UFOs. Compared to the other distributions, the UFO dis- 
tribution is much more skewed toward lower values of q with 
the bulk of UFOs having q « 0.4. The other OFGs tend 
to have a fairly large number of low q objects as well, but 
also have more q > 0.5 objects similar to what is seen in the 
full-parent sample. 


Table 3 (continued) 


Parameter UFOs OFGs 


NoTE— Reported quantity is the median and 
percentiles (16th and 84th) of the distribution 
for that parameter. 


5.2. 8D shapes from axis ratio distribution modelling 


Here we explore the intrinsic shapes of the UFOs by 
modeling the distribution of their projected axis ratios 
(e.g. Chang et al. 2013; van der Wel et al. 2014). In gen- 
eral, the shape of a galaxy can be approximated as an el- 
lipsoid characterized by three axis lengths, A > B = C, 
which can be used to define the ellipticity (E) and tri- 
axiality (T) parameters. Within the family of ellipsoids, 
there are three broad categories: the oblate ellipsoid 
(A ~ B > C, disky), prolate ellipsoid (A > B ~ C, 
elongated), or the spheroid (A ~ B ~ C). The pro- 
jected axis-ratios of these shapes can be inferred from 
how how each ellipsoid type would appear in projection 
from a set of randomly drawn viewing angles with dis- 
tinct shapes having unique projected axis-ratio distribu- 
tions (e.g. Chang et al. 2013; van der Wel et al. 2014) 
characterized by values of E and T and thus values of 
A, B, and C. For example, a q distribution peaking be- 
tween 0.8 and 1.0 would be predominately comprised 
of spheroids as they can only have large observed axis- 
ratios. A disk population would have a roughly flat q dis- 
tribution reflecting that disks can be observed edge-on, 
face-on, and everything in-between, while only observ- 
ing inclined disks would result in an axis-ratio distribu- 
tion without any higher q objects. Prolate populations 
have axis-ratio distributions that peak at shorter values 
(q ~ 0.4) with a drop at q ~ 0.2 and a tail toward higher 
axis-ratios (see left panel of Figure 11). 

In Figure 10, we plot the observed axis-ratio distri- 
butions of the full F444W-selected parent sample, all 
optically faint galaxies, and the size-selected subset (the 
UFOs). We find that the UFOs typically have q « 0.5 
and the smaller OFGs typically have q > 0.5. Taken 
at face value, this suggests that the UFOs are largely 
prolate or disk shaped and the more compact OFGs are 


more spheroid-dominated. We additionally perform a 
KS test on the UFO and OFG axis-ratio distributions to 
constrain the probability that these samples were drawn 
from the same underlying population. We find a p-value 
of ~ 1078 showing that, at the very least, the UFOs and 
OFGs represent different populations in terms of their 
axis-ratios, though we note that owing to their sizes be- 
ing closer to the scale of the PSF, the shapes of the com- 
pact OFGs are much less certain. Figure 11 shows the 
observed UFO axis ratio distribution compared to vari- 
ous mock populations. Our sample of UFOs has many 
more objects with q ~ 0.1 — 0.2 than the sample of Nel- 
son et al. (2023) and as such neither inclined thick-disk 
or prolate populations by themselves adequately match 
the observed UFO q distribution. The tail of very low 
axis-ratios can only be achieved by having a popula- 
tion of intrinsically inclined thin-disks capable of reach- 
ing such low q values. However, as can be seen in the 
right panel of Figure 11, an inclined thin-disk population 
by itself over-predicts these low q objects and under- 
predicts those with higher q. For this reason, we show 
combinations of prolate and disk-like populations on the 
right. With the exception of the prolate+thin-disk (no 
inclination restriction) which overpredicts higher q ob- 
jects, the various combinations all seem to fairly well 
match the observed distribution making it challenging 
to uniquely assign an intrinsic shape to the UFOs. 

For these reasons, we turn to a modeling scheme that 
directly fits the observed axis-ratio distributions given 
some assumptions about the true underlying distribu- 
tion. We use the BEAST (Price et al. in prep.) axis- 
ratio modeling code to perform nested sampling with 
DYNESTY (Speagle 2020) to infer the intrinsic shapes of 
the UFOs. Here, we briefly summarize the functional- 
ity of this software. The ellipticity (E = 1-C?) and 
triaxiality (T = [1-B?] / [1-C?]) can be used to define 
specific sub-types of ellipsoids characterized by ranges of 
values of E and T reflecting the underlying relationships 
between A, B, and C. In practice, one can combine an 
observed q distribution and, assuming the galaxies come 
from a population that can be described by a distribu- 
tion of E and T, perform MCMC or nested sampling of 
the posterior to determine the most likely values for the 
parameters of the E and T distributions, and thus on 
the intrinsic shapes (i.e., oblate/disky or prolate) given 
these assumptions. 

'The BEAST modeling assumes all galaxies in the sam- 
ple are drawn from a single population with Gaussian 
distributions of E and T, characterized by four free pa- 
rameters: centers yp and ur and sigmas, og and or. 
In our default model, we consider galaxies that are ran- 
domly observed from any viewing angle. However, we 
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also employ a secondary model that restricts the inclina- 
tion to be between 50—90°, and a third that additionally 
enforces oblateness (i.e., diskiness) by restricting the tri- 
axiality parameter « 0.33, in addition to restricting the 
inclination. These secondary models encode the physical 
expectation that these objects are inclined disks. Each 
model fits for the q distribution in 10 bins. For each 
parameter, we assume uniform priors from zero to unity 
except for the third model where we restrict the range 
of T to be between zero and 0.33. In what follows, we 
define four unique regions of intrinsic axis-ratio space (C 
versus B) that represent distinct shapes and can be seen 
in the middle panel of Figure 12. These are spheroidal 
(green/upper right), prolate (blue/lower left), thin-disk 
(red/bottom of circular region), and thick-disk (pur- 
ple/top of circular region). The results of fitting these 
three models are summarized in Figure 12 where we 
show the best-fit (maximum a posteriori) model q dis- 
tribution, random draws from the distribution of the 
best-fit model in intrinsic axis-ratio space, and the frac- 
tions of each shape that can be attributed to that single 
population model. 

From the top panels in Figure 12, we can see that each 
different version of the model (free, inclined, inclined 
and oblate) does a reasonably good job at reproduc- 
ing the observed UFO axis-ratio distribution. However, 
the model with no triaxiality or inclination constraints 
achieves this agreement in a markedly different way than 
the constrained models. In particular, when there are 
no constraints on the inclination or triaxiality of the 
population, the model prefers a population that is dom- 
inated by intrinsically prolate objects (~ 90%) with only 
a ~ 10% contribution from disks (mostly thin) and es- 
sentially zero spheroidal objects. This tells us that the 
UFO q distribution can be generated from a popula- 
tion that is dominated by prolate ellipsoids with a small 
contribution from thin-disks to reproduce the lowest q 
UFOs. 

In the other two models, the story is very different. 
These disk assumption models can also produce axis- 
ratio distributions that are fairly consistent with the 
data showing that disk-based distributions are not in- 
consistent with the observations. In the middle panel 
of the bottom row, we can see that the restricted incli- 
nation model still predicts that prolate objects are the 
dominant contributor to the observed q distribution, but 
only slightly with a ~ 40% contribution from thick and 
thin disks. Thus, even if we assume that we are only 
observing galaxies closer to edge-on, the axis-ratio mod- 
eling still statistically prefers a contribution from prolate 
objects. Moving to the third panel of the bottom row, 
we can see that if we restrict the triaxiality to be less 
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Figure 11. Left panel: Observed UFO axis-ratio distribution (black) compared with different q distributions generated for 
different mock galaxy populations (prolate/blue, thick-disk/solid orange, thin-disk/dashed orange) showing the plausible galaxy 
types that could give rise to the observed quro distribution. Right panel: Same as the left panel, but for different linear 
combinations of distinct populations (inclined thick/thin disks in orange, prolate/thin-disk in blue, and prolate/inclined thin 
disk in purple). It is difficult to reproduce quro without both an inclined thin-disk population and a prolate population. 


than 0.33 as well as restricting the inclination, then the 
model predicts the dominant contribution coming from 
thick and thin disks (~ 60%) with a contribution from 
prolate ellipsoids of at most ~ 40%. This shows that 
even in the most disk-friendly modeling scenario, there 
is evidence for a contribution from prolate set of objects. 
That some optically faint objects might be prolate is a 
surprising finding as it is not straightforward how a pro- 
late object viewed along its long-axis could have dust 
column densities large enough to render the galaxy op- 
tically faint. We can restrict the triaxiality even more, 
which will put more objects into the disk category, but 
this leads to much worse agreement between observed 
and model g distributions and requires a ~ 20% spheroid 
contribution. 

We note here a few modeling caveats that should be 
mentioned before we continue to discuss broader impli- 
cations of these results. First, we assume that the in- 
trinsic shapes of the UFOs can be accurately modeled 
as ellipsoids, which is a good approximation, but will 
clearly break down for certain objects. For example, 
ongoing galaxy mergers, galaxies with clumpy features, 
or galaxies with spiral arms and/or dust lanes could de- 
viate somewhat from the ellipsoidal shape assumption. 
In fact, in some of our UFOs (Figure 2), we see evidence 
of clumps (29312), spiral features (187160), or multi- 
ple components / possible mergers (e.g., 167032 and 
172813). These non-smooth components could impact 


the ensemble axis-ratio distribution leading to higher or 
lower inferences on the fractions of one or more shape 
subcategories. Additionally, if the UFOs are not a uni- 
fied population of intrinsic shapes that can be character- 
ized by a single distribution, then our modeling assump- 
tions start to break down and the inferred parameters 
are not necessarily telling us about the true breakdown 
of intrinsic shapes. 


6. DISCUSSION 


This paper presents an analysis of the stellar popu- 
lations and structural properties of an unexpected pop- 
ulation of elongated, optically faint galaxies from the 
JADES survey. Here, we discuss why these results are 
surprising in terms of existing studies of the intrinsic 
shapes of galaxies, what drives the optical faintness of 
these objects, and their subsequent evolution. 

Probably the most surprising aspect of this popula- 
tion of objects is their structure: they have fairly large 
radii, nearly exponential light distributions, and are 
highly elongated. Beginning with their sizes, the left- 
most panel of Figure 9 shows that a population of op- 
tically faint galaxies exist (which we call UFOs) whose 
sizes are larger than average star-forming galaxies at 
their masses and redshifts. This is surprising because at 
a fixed dust mass, galaxies with a dust distribution that 
is more concentrated will have higher dust column den- 
sities and hence larger values of dust attenuation (e.g. 
Nelson et al. 2014). These higher dust column densities 
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Figure 12. Summary of the results from our axis-ratio modeling with the BEAST. The top row shows the observed axis-ratio 
distributions in black compared with the best-fitting (maximum a posteriori) model distributions in red. Lighter red lines 
show q distributions for random draws of the posterior distributions of ug, ur,cE,or. The middle row shows random draws 
from the best-fitting model distribution in intrinsic axis-ratio space illustrating how the four distinct regions are defined and 
what the contribution from each region is. In the bottom row, we show the predicted fractions of four different galaxy shapes: 
elongated/prolate (blue), spheroidal (green), thick disk (purple) and thin disk (red). From left to right, we show the fully free 
version of the axis-ratio model, the model with restricted inclination, and the model with restricted inclination and a stronger 
preference for oblate galaxies. Here, vertical lines are the best-fit value for the fraction, horizontal lines show the 68% confidence 
interval, and the violin plots (shaded areas) give a representation of the distribution with wider regions having more probability. 
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could come from compact sizes or concentrated distri- 
butions (i.e., high Sérsic indices). As such, we may ex- 
pect optically faint galaxies to preferentially have small 
sizes and / or high Sérsic indices. Thus, while the naive 
physical expectation for the structures of optically faint 
galaxies would be to have compact sizes to drive their 
high values of dust attenuation, this does not appear to 
be the case with the distributions of the optically-faint 
and bluer F444W selected parent samples being simi- 
lar. Specifically, amongst objects with stellar masses of 
greater than 101? M, sizes > 2kpc, 39% are optically 
faint while amongst objects with sizes less than 2kpc, 
2296 are optically faint. Optical faintness does not ap- 
pear to correlate with size. This is surprising given the 
findings of Gómez-Guijarro et al. (2023) who find that 
amongst star forming galaxies at 3 < z « 7.5, those with 
Ay > 1 have ~ 30% smaller sizes. Although, recent sim- 
ulations find that orientation (i.e., viewing angle) and 
not physical characteristics, is the main determinate of 
whether a galaxy will be optically faint or not, with 
galaxies only appearing optically-faint in some fraction 
(depending on the object) of the orientations (Cochrane 
et al. 2023). 

The second surprise is their Sérsic indices: the bulk of 
the UFOs and the optically faint galaxies in general have 
n « 2 (see Fig. 9). Previous multi-wavelength studies 
looking at the morphologies of dusty galaxies in both 
optical/UV and sub-mm wavelengths general find that 
these sources are more compact at longer observed wave- 
lengths (e.g., Hodge et al. 2016; Tadaki et al. 2017; Nel- 
son et al. 2019) suggesting that these galaxies may be in 
the process of building dense stellar bulges (e.g., Lang 
et al. 2014; Whitaker et al. 2017) with dust attenua- 
tion causing these bulges to be missed at shorter wave- 
lengths. Bulge-dominated galaxies are typically char- 
acterized by high Sérsic indices reflecting more centrally 
concentrated light profiles, which is the opposite of what 
we see in the sample of optically faint galaxies studied in 
this paper. The UFOs and optically-faint galaxies in our 
sample have light distributions more indicative of expo- 
nential profiles suggesting they are strongly disk domi- 
nated (though prolate is also a possibility). The massive 
portion of the galaxies in our sample are approaching 
stellar masses at which they will likely quench yet most 
show little evidence for having or building bulges which 
appears to be a prerequisite for quenching (e.g. Lang 
et al. 2014; Whitaker et al. 2017). Although most of the 
UFOs have low Sérsic indices in F444W emission, it is 
possible that even the F444W light is attenuated in the 
central regions and a bulge component could appear in 
even redder wavelengths. MIRI or ALMA data would 


be needed in order to definitively rule out the presence 
of a highly obscured bulge component. 

'The third surprise is that most of the UFOs are elon- 
gated. These objects are selected based on their fluxes, 
colors, and sizes; there is no explicit selection for them 
to be elongated. The most obvious explanation per- 
haps for a highly reddened cohort of large galaxies to 
be elongated is that we are looking at an edge-on subset 
of a disk galaxy population. Another possible explana- 
tion for high dust attenuation is geometric: if galaxies 
are intrinsically disky, then the dust column density will 
be much larger when viewed edge-on than face-on (e.g. 
Wild et al. 2011; Patel et al. 2012). If the UFOs are 
in fact disk-dominated objects, viewing them edge-on 
could explain why they are optically faint and would 
not have appeared in previous studies. 

While there is significant physical reason to believe 
the UFOs are intrinsically edge-on disks, axis-ratio mod- 
elling suggests a significant fraction of objects may be 
intrinsically prolate instead of oblate/disky. This type 
of statistical modeling was not possible in Nelson et al. 
(2023) given the small sample size; having the large sam- 
ple provided by JADES is essential. In a framework 
assuming UFOs are a single population, axis-ratio mod- 
elling suggests that at least 3096 of UFOs are intrinsi- 
cally prolate. This is very surprising. In an intrinsi- 
cally prolate population, a physical argument akin to 
the higher dust columns in inclined disk galaxies, would 
suggest a preference for optical faintness when looking 
at the prolate object with the long axis along the line 
of sight. In this configuration, the object would ap- 
pear small and circular (perhaps like some of the smaller 
optically-faint galaxies in our sample) and it would phys- 
ically make sense for these objects to have large dust 
columns. Viewed with the long axis perpendicular to 
the line of sight, it is harder to understand how dust 
attenuation could be strong enough to produce optical 
faintness. 

'To put these results in context, we compare to previ- 
ous studies of the intrinsic shapes of galaxies as a func- 
tion of redshift and stellar mass. Our stellar popula- 
tion modelling suggests that these objects have fairly 
high stellar masses; observations of star-forming galax- 
ies (SFGs) with HST over the redshift range 1 < z < 3 
have shown that high mass (M, > 10!? M5) SFGs are 
predominately disks (e.g., van der Wel et al. 2014; Zhang 
et al. 2019) especially at z « 2. However, at higher 
redshifts (z = 2) and lower masses (M, < 109? M) 
the shapes of SFGs become dominated by prolate ob- 
jects. In other words, the fraction of star-forming pro- 
late galaxies increases towards higher redshifts and lower 
masses, while the opposite is true of star-forming disks. 


mmm All Galaxies (M « > 1019 ^M,) 
All Galxies (M« < 1019 4M,) 
[==] UFOs (M. > 101?^M,) 
"I! UFOs (M. < 101%4Mo) 


Figure 13. Axis-ratio distributions for the UFOs and the 
entire sample of galaxies considered in this work split into 
low mass (M, < 10'°*Mo) and high mass (M, < 10'°*Mo) 
samples. It can clearly be seen that the majority of the low q 
objects in our whole sample are predominately coming from 
lower mass UFOs, but with a strong contribution from the 
higher mass UFOs as well. The q distribution of low-mass 
UFOs is largely consistent with a prolate population, but the 
lower sample size of the high-mass group make it harder to 
come to strong conclusions about the intrinsic shapes of the 
higher mass UFOs. 


For this reason, and given the fairly wide range of in- 
ferred UFO stellar masses, we consider the possibility 
that the UFOs are composed of a lower mass popula- 
tion and a higher mass population with distinct axis- 
ratio distributions. Figure 13 shows the axis ratio dis- 
tributions of UFOs and the parent population of F444W 
selected galaxies divided by mass at M, = 10!°4Mo. 
As expected, the axis-ratio distribution of the low mass 
UFOs is more skewed towards lower q values than that 
of the high mass UFOs and axis-ratio modeling suggests 
nearly 100% are prolate. There are not enough galaxies 
in the high mass bin to determine their intrinsic shapes 
from statistical modelling but their axis-ratio distribu- 
tion is certainly not that of a purely prolate population. 

It has been shown in Vega-Ferrero et al. (2023), based 
on a machine-learning based classification scheme for in- 
ferring galaxy morphologies, that approximately half of 
galaxies visual classified as disks from CEERS imaging 
are more consistent with prolate or spheroidal popu- 
lations than being true disks. Thus, we may indeed 
be looking at some objects that visually appear to be 
disks, but in reality are prolate. Further, it is worth 
noting that populations of massive, but quenched pro- 
late galaxies at z ~ 0 have been identified in cosmo- 
logical simulations (e.g., Ebrová & Lokas 2017; Thob 
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et al. 2019) as well as observationally (e.g., Tsatsi et al. 
2017). Ebrová & Lokas (2017) suggests that these pro- 
late objects are either formed prolate (or become pro- 
late early in their lifetimes) or formed through radially 
aligned mergers that induce the stellar shapes towards 
prolate ellipsoids. These types of observations could 
be consistent with our plausible finding of star-forming 
prolate intermediate mass (M, < 10!°4M, galaxies at 
2 < z < 6 which could then be the progenitors of mas- 
sive quenched galaxies at z ~ 0. Although we can rule 
out that the UFOs are not a bulge-dominated popula- 
tion in F444W light, definitively determining the intrin- 
sic shapes of these objects will require kinematic mea- 
surements. 

'Turning back to why these objects are optically-faint, 
in Figure 6 we show the distribution of optically faint 
galaxies and our F444W-selected parent sample in a 
plane of two likely culprits: age and dust attenuation, 
which both act to redden the SEDs of galaxies. The 
most significant effect is dust attenuation: most galaxies 
at 2< z < 6 and F444W < 24.5 that have A, > 2 are 
optically faint. There is a secondary trend with age in 
which the optically faint galaxies selected from JADES 
are on average slightly older than the parent population 
but this is decidedly subdominant to the dust attenua- 
tion. Further, as shown in Fig. 9, 58% of galaxies in our 
F444W-selected parent sample with q « 0.3 are optically 
faint. On the other hand, both size and Sérsic index 
do not appear driving optical faintness (at least as de- 
fined in this paper). Thus, having a low projected axis- 
ratio appears to be the best predictor of optical faintness 
amongst the standard morphological indicators — effec- 
tive radius, Sérsic index, and axis-ratio. That said, a 
very small fraction of the compact optically-faint galax- 
ies (i.e. non-UFOs) have q « 0.3 so inclination does not 
appear to provide an explanation for the optical faint- 
ness of these more compact objects. 

Finally, we consider how this enigmatic galaxy popu- 
lation may evolve toward lower redshifts. With stellar 
masses of log(M..)> 10 at z ~ 2 — 3, UFOs are likely 
to evolve into quiescent galaxies with log(M,)> 11 in 
the local universe (e.g. Behroozi et al. 2019). In fact, in 
Figure 5, we can see that a small fraction of our three 
samples have SFRs up to ^ 3 dex below the z — 2.5 
SFMS from Leja et al. (2022) suggesting that these 
galaxies have already ceased their star-formation. This 
is a sharp transition to quiescence at a characteristic 
mass of log(M..)~ 10.4 consistent with previous studies 
finding that quenching is efficient at a similar mass (e.g., 
Contini et al. 2020). With the typically small quantity 
of dust attenuation in quiescent galaxies in the local 
universe, the dust attenuation seen in these objects at 
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z > 2 will need to almost completely disappear in the 
intervening time (e.g. Whitaker et al. 2021). Significant 
structural change is also likely as they need to transition 
from prolate or oblate morphologies to having significant 
to dominant bulge components. 


7. CONCLUSIONS 


We have performed a morphological and stellar popu- 
lation analysis of a sample of 112 optically-faint galaxies 
with a specific focus on the elongated sub-set of optically 
faint galaxies, dubbed UFOs. Specifically, we explored 
how these galaxies are situated in certain fundamental 
galaxy scaling relations such as the size-mass plane and 
the star-forming main-sequence as well as investigations 
into how the stellar populations vary between the inner 
and outer regions and finally an in depth look at the 
morphological parameters of these objects such as their 
sizes, Sérsic indices, and axis-ratios. Our main goal has 
been to place stronger constraints on the nature of the 
UFO galaxy population and to place the optically faint 
and UFO populations within the larger context of the 
massive, star-forming galaxy population. Our main re- 
sults are summarized below. 


1. We identify 112 optically-faint objects in the 
JADES survey of which 56 are classified as UFOs 
with apparent sizes larger than 0.25". The UFOs 
are typically between 2 « z « 4 with a range 
of stellar masses (M, ~ 109?-!! M), high SFRs, 
and high Ay. Compared to a mass and redshift 
matched sample of bluer F444W-selected galax- 
ies in JADES, the UFOs tend to have higher stel- 
lar masses and SFRs, but have much lower stellar 
masses than the Herschel detected ultrared galax- 
ies from Ma et al. (2019). 


2. UFOs have red colors throughout the extent of 
their bodies, but still possess slight negative color 
gradients with outskirts that are less red than the 
interiors. These color gradients are likely driven 
by increased central dust concentration and not 
by the ages of the stellar populations. 


3. The observed UFO  axis-ratio distribution in 
F444W is consistent with observing some combi- 
nation of randomly oriented disks and randomly 
oriented prolate galaxies. We perform a detailed 
Bayesian modeling of the q distribution using 
the BEAST and find that the UFOs can not be 
uniquely said to consist of one single galaxy shape. 
It is possible that the UFOs consist of a lower mass 
prolate population and a higher mass oblate/disk 
population, but without a larger sample it is hard 
to verify this. 


4. We find that the strongest predictors of optical 
faintness in the galaxy populations are Ay and q. 
That increasing Ay is a strong predictor of op- 
tical faintness is not surprising given that larger 
dust attenuation will naturally lead to less opti- 
cal light escaping a galaxy. Morphologically, it is 
somewhat surprising that the strongest predictor 
for optical faintness is a low q as it was thought 
that small physical size would be the easiest way to 
drive large enough dust columns to make galaxies 
optically-faint, but that is not what we see given 
the large sizes of UFOs. 


UFOs are an unanticipated galaxy population that 
JWST has made easy to reveal and study. The exis- 
tence of optically-faint galaxies with extended structure 
was surprising owing to the expectation that only com- 
pact sizes could lead to large enough dust columns and 
as early JWST identified optically-faint galaxies were 
found to be relatively compact (Gómez-Guijarro et al. 
2023). The suggestion that a significant fraction of these 
objects are prolate instead of disks is an interesting pos- 
sibility with implications for the z > 2 shapes of star- 
forming galaxies and how these shapes evolve to the 
present day. Future work will increase these samples 
to larger numbers and seek stronger constraints on the 
morphologies with kinematics from JWST spectra al- 
lowing us to come to a definitive conclusion on the in- 
trinsic 3D shapes of UFOs. 


J.L.G gratefully acknowledges support provided by 
NASA through grants 20-ASTRO20-0200 and HST- 
AR-16146 and support from the NRAO through grant 
SOSPADA-025. E.J.N acknowledges support provided 
by NASA through grants HST-AR-16146 and JWST- 
GO-01895. The research of C.C.W is supported by 
NOIRLab, which is managed by the Association of 
Universities for Research in Astronomy (AURA) un- 
der a cooperative agreement with the National Science 
Foundation. S.H.P acknowledges support through grant 
JWST-GO-2561. The Cosmic Dawn Center (DAWN) is 
funded by the Danish National Research Foundation un- 
der grant No. 140. A.J.B acknowledge funding from the 
FirstGalaxies Advanced Grant from the European Re- 
search Council (ERC) under the European Union's Hori- 
zon 2020 research and innovation programme (Grant 
agreement No. 789056). W.B. acknowledges support by 
the Science and Technology Facilities Council (STFC), 
ERC Advanced Grant 695671 "QUENCH". K.B. ac- 
knowledges support from the Australian Research Coun- 
cil Centre of Excellence for All Sky Astrophysics in 


3 Dimensions (ASTRO 3D), through project number 
CE170100013. E.C.L acknowledges support of an STFC 
Webb Fellowship (ST/W001438/1). D.J.E is supported 
as a Simons Investigator and by JWST/NIRCam con- 
tract to the University of Arizona, NAS5-02015. R.H 
acknowledges support from the Johns Hopkins Univer- 
sity, Institute for Data Intensive Engineering and Sci- 
ence (IDIES). R.M. acknowledges support by the Sci- 
ence and Technology Facilities Council (STFC), by the 
ERC through Advanced Grant 695671 QUENCH, and 
by the UKRI Frontier Research grant RISEandFALL. 
R.M. also acknowledges funding from a research profes- 
sorship from the Royal Society. B.D.J., G.R., M.R., and 
B.E.R. acknowledge support from a JWST/NIRCam 


21 


contract to the University of Arizona NAS5-02015. This 
work utilized the Alpine high performance computing 
resource at the University of Colorado Boulder. Alpine 
is jointly funded by the University of Colorado Boul- 
der, the University of Colorado Anschutz, and Colorado 
State University 


Facilities: JWST(NIRCam) 


Software:  astropy (Astropy Collaboration et al. 
2013, 2018, 2022), DvNESTY (Speagle 2020), GAL- 
FIT (Peng et al. 2002, 2010), LENSTRONOMY (Birrer & 
Amara 2018), matplotlib (Hunter 2007), numpy (Harris 
et al. 2020), photutils (Bradley et al. 2020), PRosPEC- 
TOR (Johnson et al. 2021) 


REFERENCES 


Anderson, J. 2016, Empirical Models for the WFC3/IR 
PSF, Instrument Science Report WFC3 2016-12, 42 pages 
Anderson, J., & King, I. R. 2000, PASP, 112, 1360, 
doi: 10.1086/316632 
Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., 
et al. 2013, A&A, 558, A33, 
doi: 10.1051/0004-6361/201322068 
Astropy Collaboration, Price-Whelan, A. M., Sipócz, B. M., 
et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f 
Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., 
et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74 
Barger, A. J., Cowie, L. L., Blair, A. H., & Jones, L. H. 
2022, ApJ, 934, 56, doi: 10.3847/1538-4357/ac67e7 
Barger, A. J., Cowie, L. L., Sanders, D. B., et al. 1998, 
Nature, 394, 248, doi: 10.1038/28338 
Barger, A. J., Wang, W. H., Cowie, L. L., et al. 2012, ApJ, 
761, 89, doi: 10.1088/0004-637X /761/2/89 
Barro, G., Kriek, M., Pérez-González, P. G., et al. 2016, 
ApJL, 827, L32, doi: 10.3847/2041-8205/827/2/L32 
Barrufet, L., Oesch, P. A., Weibel, A., et al. 2023, MNRAS, 
doi: 10.1093/mnras/stad947 
Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 
2019, MNRAS, 488, 3143, doi: 10.1093/mnras/stz1182 
Bell, E. F., van der Wel, A., Papovich, C., et al. 2012, ApJ, 
753, 167, doi: 10.1088/0004-637X/753/2/167 
Birrer, S., & Amara, A. 2018, Physics of the Dark Universe, 
22, 189, doi: 10.1016/j.dark.2018.11.002 
Birrer, S., Shajib, A., Gilman, D., et al. 2021, The Journal 
of Open Source Software, 6, 3283, 
doi: 10.21105/joss.03283 
Bradley, L., Sipócz, B., Robitaille, T., et al. 2020, 
astropy/photutils: 1.0.1, 1.0.1, Zenodo, Zenodo, 
doi: 10.5281/zenodo.4049061 


Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, 
ApJ, 686, 1503, doi: 10.1086/591786 

Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2022, 
JWST Calibration Pipeline, 1.8.1, Zenodo, 
doi: 10.5281 /zenodo.7215690 

Byler, N., Dalcanton, J. J., Conroy, C., & Johnson, B. D. 
2017, ApJ, 840, 44, doi: 10.3847/1538-4357/aa6c66 

Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 
533, 682, doi: 10.1086/308692 

Casey, C. M., Narayanan, D., & Cooray, A. 2014, PhR, 
541, 45, doi: 10.1016/j.physrep.2014.02.009 

Chabrier, G. 2003, PASP, 115, 763, doi: 10.1086/376392 

Chan, J. C. C., Beifiori, A., Mendel, J. T., et al. 2016, 
MNRAS, 458, 3181, doi: 10.1093/mnras/stw502 

Chang, Y.-Y., van der Wel, A., Rix, H.-W., et al. 2013, 
ApJ, 773, 149, doi: 10.1088 /0004-637X/773/2/149 

Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718, 
doi: 10.1086/309250 

Chen, C.-C., Cowie, L. L., Barger, A. J., Wang, W.-H., & 
Williams, J. P. 2014, ApJ, 789, 12, 
doi: 10.1088/0004-637X/789/1/12 

Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102, 
doi: 10.3847 /0004-637X /823/2/102 

Cochrane, R. K., Anglés-Alcázar, D., Cullen, F., & 
Hayward, C. C. 2023, arXiv e-prints, arXiv:2310.08829, 
doi: 10.48550/arXiv.2310.08829 

Conroy, C., Gunn, J. E., & White, M. 2009, ApJ, 699, 486, 
doi: 10.1088/0004-637X /699/1/486 

Conroy, C., White, M., & Gunn, J. E. 2010, ApJ, 708, 58, 
doi: 10.1088/0004-637X /708/1/58 

Contini, E., Gu, Q., Ge, X., et al. 2020, ApJ, 889, 156, 
doi: 10.3847/1538-4357 /ab6730 


22 


Coppin, K., Chapin, E. L., Mortier, A. M. J., et al. 2006, 
MNRAS, 372, 1621, 
doi: 10.1111/j.1365-2966.2006.10961.x 

Cowie, L. L., González-López, J., Barger, A. J., et al. 2018, 
ApJ, 865, 106, doi: 10.3847/1538-4357 /aadc63 

Draine, B. T., & Li, A. 2007, ApJ, 657, 810, 
doi: 10.1086/511055 

Ebrová, I., & Lokas, E. L. 2017, ApJ, 850, 144, 
doi: 10.3847/1538-4357 /aa96ff 

Eisenstein, D. J., Willott, C., Alberts, S., et al. 2023, arXiv 
e-prints, arXiv:2306.02465, 
doi: 10.48550/arXiv.2306.02465 

Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 
533, A119, doi: 10.1051/0004-6361/201117239 

Falcón-Barroso, J., Sánchez-Blázquez, P., Vazdekis, A., 
et al. 2011, A&A, 532, A95, 
doi: 10.1051/0004-6361/201116842 

Fasano, G., & Franceschini, A. 1987, MNRAS, 225, 155, 
doi: 10.1093/mnras/225.1.155 

Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, 
PASP, 110, 761, doi: 10.1086/316190 

Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, 
RMxAA, 53, 385, doi: 10.48550/arXiv.1705.10877 

Franco, M., Elbaz, D., Béthermin, M., et al. 2018, A&A, 
620, A152, doi: 10.1051/0004-6361/201832928 

Franx, M., van Dokkum, P. G., Fórster Schreiber, N. M., 
et al. 2008, ApJ, 688, 770, doi: 10.1086/592431 

Fujimoto, S., Ouchi, M., Shibuya, T., & Nagai, H. 2017, 
ApJ, 850, 83, doi: 10.3847/1538-4357/aa93e6 

Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, 
PASP, 135, 068001, doi: 10.1088/1538-3873/acd1b5 

Gómez-Guijarro, C., Elbaz, D., Xiao, M., et al. 2022, A&A, 
658, A43, doi: 10.1051/0004-6361/202141615 

Gómez-Guijarro, C., Magnelli, B., Elbaz, D., et al. 2023, 
arXiv e-prints, arXiv:2304.08517, 
doi: 10.48550/arXiv.2304.08517 

Gottumukkala, R., Barrufet, L., Oesch, P. A., et al. 2023, 
arXiv e-prints, arXiv:2310.03787, 
doi: 10.48550/arXiv.2310.03787 

Gullberg, B., Smail, I., Swinbank, A. M., et al. 2019, 
MNRAS, 490, 4956, doi: 10.1093/mnras/stz2835 

Guo, Y., Giavalisco, M., Cassata, P., et al. 2011, ApJ, 735, 
18, doi: 10.1088/0004-637X /735/1/18 

Hainline, K. N., Johnson, B. D., Robertson, B., et al. 2023, 
arXiv e-prints, arXiv:2306.02468, 
doi: 10.48550/arXiv.2306.02468 

Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 
2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2 

Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 
208, 19, doi: 10.1088/0067-0049/208/2/19 


Hodge, J. A., Swinbank, A. M., Simpson, J. M., et al. 2016, 
ApJ, 833, 103, doi: 10.3847/1538-4357/833/1/103 

Hsu, L.-Y., Cowie, L. L., Barger, A. J., & Wang, W.-H. 
2017, ApJ, 850, 189, doi: 10.3847/1538-4357/aa9619 

Hughes, D. H., Serjeant, S., Dunlop, J., et al. 1998, Nature, 
394, 241, doi: 10.1038/28328 

Hunter, J. D. 2007, Computing in Science and Engineering, 
9, 90, doi: 10.1109/M CSE.2007.55 

Illingworth, G., Magee, D., Bouwens, R., et al. 2016, arXiv 
e-prints, arXiv:1606.00841, 
doi: 10.48550/arXiv.1606.00841 

Ji, Z., Williams, C. C., Tacchella, S., et al. 2023, arXiv 
e-prints, arXiv:2305.18518, 
doi: 10.48550/arXiv.2305.18518 

Johnson, B. D., Leja, J., Conroy, C., & Speagle, J. S. 2021, 
ApJS, 254, 22, doi: 10.3847/1538-4365 /abef67 

Kokorev, V., Jin, S., Magdis, G. E., et al. 2023, ApJL, 945, 
L25, doi: 10.3847/2041-8213/acbd9d 

Lang, P., Wuyts, S., Somerville, R. S., et al. 2014, ApJ, 
788, 11, doi: 10.1088/0004-637X /788/1/11 

Leja, J., Carnall, A. C., Johnson, B. D., Conroy, C., & 
Speagle, J. S. 2019a, ApJ, 876, 3, 
doi: 10.3847/1538-4357 /ab133c 

Leja, J., Johnson, B. D., Conroy, C., & van Dokkum, P. 
2018, ApJ, 854, 62, doi: 10.3847/1538-4357/aaa8db 

Leja, J., Johnson, B. D., Conroy, C., van Dokkum, P. G., & 
Byler, N. 2017, ApJ, 837, 170, 
doi: 10.3847/1538-4357 /aabffe 

Leja, J., Johnson, B. D., Conroy, C., et al. 2019b, ApJ, 877, 
140, doi: 10.3847/1538-4357/abld5a 

Leja, J., Speagle, J. S., Ting, Y.-S., et al. 2022, ApJ, 936, 
165, doi: 10.3847/1538-4357/ac887d 

Ma, J., Cooray, A., Nayyeri, H., et al. 2019, ApJS, 244, 30, 

doi: 10.3847/1538-4365/ab4194 

Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415, 

doi: 10.1146/annurev-astro-081811-125615 

Manning, S. M., Casey, C. M., Zavala, J. A., et al. 2022, 

ApJ, 925, 23, doi: 10.3847/1538-4357 /ac366a. 

Marrone, D. P., Spilker, J. S., Hayward, C. C., et al. 2018, 

Nature, 553, 51, doi: 10.1038/nature24629 

Mathews, E. P., Leja, J., Speagle, J. S., et al. 2023, arXiv 
e-prints, arXiv:2306.16442, 
doi: 10.48550/arXiv.2306.16442 

Miller, T. B., van Dokkum, P., & Mowla, L. 2023, ApJ, 
945, 155, doi: 10.3847/1538-4357 /acbc74 

Mosleh, M., Tacchella, S., Renzini, A., et al. 2017, ApJ, 
837, 2, doi: 10.3847/1538-4357 /aabf14 

Nelson, E., van Dokkum, P., Franx, M., et al. 2014, Nature, 
513, 394, doi: 10.1038/nature13616 


Nelson, E. J., van Dokkum, P. G., Fórster Schreiber, N. M., 

et al. 2016, ApJ, 828, 27, 

doi: 10.3847/0004-637X /828/1/27 

Nelson, E. J., Tadaki, K.-i., Tacconi, L. J., et al. 2019, ApJ, 

870, 130, doi: 10.3847/1538-4357 /aaf38a 

Nelson, E. J., Suess, K. A., Bezanson, R., et al. 2023, 

ApJL, 948, L18, doi: 10.3847/2041-8213/acclel 

Nenkova, M., Sirocky, M. M., Ivezié, Z., & Elitzur, M. 

2008a, ApJ, 685, 147, doi: 10.1086/590482 

Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezić, Z., & 

Elitzur, M. 2008b, ApJ, 685, 160, doi: 10.1086/590483 

Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 
507, 1793, doi: 10.1051/0004-6361/200912497 

Oke, J. B. 1974, ApJS, 27, 21, doi: 10.1086/190287 

Oteo, L, Ivison, R. J., Dunne, L., et al. 2016, ApJ, 827, 34, 
doi: 10.3847/0004-637X /827 /1/34 

Patel, S. G., Holden, B. P., Kelson, D. D., et al. 2012, 
ApJL, 748, L27, doi: 10.1088/2041-8205/748/2/L27 

Peacock, J. A. 1983, MNRAS, 202, 615, 
doi: 10.1093/mnras/202.3.615 

Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, 
AJ, 124, 266, doi: 10.1086/340952 

—. 2010, AJ, 139, 2097, doi: 10.1088/0004-6256/139/6/2097 

Pérez-González, P. G., Barro, G., Annunziatella, M., et al. 
2023, ApJL, 946, L16, doi: 10.3847/2041-8213/acb3a5 

Riechers, D. A., Bradford, C. M., Clements, D. L., et al. 
2013, Nature, 496, 329, doi: 10.1038/nature12050 

Rieke, M. J., Kelly, D. M., Misselt, K., et al. 2023a, PASP, 
135, 028001, doi: 10.1088/1538-3873/acac53 

Rieke, M. J., Robertson, B., Tacchella, S., et al. 2023b, 
ApJS, 269, 16, doi: 10.3847/1538-4365 /acf44d 

Rodighiero, G., Bisigello, L., Iani, E., et al. 2023, MNRAS, 
518, L19, doi: 10.1093/mnrasl/slac115 

Rujopakarn, W., Dunlop, J. S., Rieke, G. H., et al. 2016, 
ApJ, 833, 12, doi: 10.3847/0004-637X /833/1/12 

Sánchez-Blázquez, P., Peletier, R. F., Jiménez-Vicente, J., 
et al. 2006, MNRAS, 371, 703, 
doi: 10.1111/j.1365-2966.2006.10699.x 

Speagle, J. S. 2020, MNRAS, 493, 3132, 
doi: 10.1093/mnras/staa278 

Suess, K. A., Kriek, M., Price, S. H., & Barro, G. 2019a, 
ApJ, 877, 103, doi: 10.3847/1538-4357/ablbda 

—. 2019b, ApJL, 885, L22, doi: 10.3847/2041-8213/ab4db3 

—. 2021, ApJ, 915, 87, doi: 10.3847/1538-4357/abfle4 

Suess, K. A., Bezanson, R., Nelson, E. J., et al. 2022, 
ApJL, 937, L33, doi: 10.3847/2041-8213/ac8e06 

Sun, F., Egami, E., Rawle, T. D., et al. 2021, ApJ, 908, 
192, doi: 10.3847/1538-4357 /abd6e4 

Szomoru, D., Franx, M., van Dokkum, P. G., et al. 2013, 

ApJ, 763, 73, doi: 10.1088/0004-637X/763/2/73 


23 


Tacchella, S., Carollo, C. M., Forster Schreiber, N. M., 
et al. 2018, ApJ, 859, 56, doi: 10.3847/1538-4357 /aabf8b 

Tadaki, K.-i., Genzel, R., Kodama, T., et al. 2017, ApJ, 
834, 135, doi: 10.3847/1538-4357/834/2/135 

Tadaki, K.-i., Belli, S., Burkert, A., et al. 2020, ApJ, 901, 
74, doi: 10.3847/1538-4357 /abaf4a 

Thob, A. C. R., Crain, R. A., McCarthy, I. G., et al. 2019, 
MNRAS, 485, 972, doi: 10.1093/mnras/stz448 

Toft, S., Smoléié, V., Magnelli, B., et al. 2014, ApJ, 782, 
68, doi: 10.1088/0004-637X/782/2/68 

Tortora, C., Napolitano, N. R., Cardone, V. F., et al. 2010, 
MNRAS, 407, 144, doi: 10.1111/j.1365-2966.2010.16938.x 

Tsatsi, A., Lyubenova, M., van de Ven, G., et al. 2017, 
A&A, 606, A62, doi: 10.1051/0004-6361/201630218 

Umehata, H., Smail, I., Swinbank, A. M., et al. 2020, A&A, 
640, L8, doi: 10.1051/0004-6361/202038146 

van der Wel, A., Chang, Y.-Y., Bell, E. F., et al. 2014, 
ApJL, 792, L6, doi: 10.1088/2041-8205/792/1/L6 

Vega-Ferrero, J., Huertas-Company, M., Costantin, L., 
et al. 2023, arXiv e-prints, arXiv:2302.07277, 
doi: 10.48550/arXiv.2302.07277 

Walter, F., Decarli, R., Carilli, C., et al. 2012, Nature, 486, 
233, doi: 10.1038/nature11073 

Walter, F., Decarli, R., Aravena, M., et al. 2016, ApJ, 833, 
67, doi: 10.3847/1538-4357/833/1/6T 

Wang, B., Leja, J., Bezanson, R., et al. 2023, ApJL, 944, 
L58, doi: 10.3847/2041-8213/acba99 

Wang, T., Elbaz, D., Schreiber, C., et al. 2016, ApJ, 816, 
84, doi: 10.3847/0004-637X /816/2/84 

Wang, T., Schreiber, C., Elbaz, D., et al. 2019, Nature, 572, 
211, doi: 10.1038/s41586-019-1452-4 

Wang, W.-H., Barger, A. J., & Cowie, L. L. 2012, ApJ, 744, 
155, doi: 10.1088/0004-637X/744/2/155 

Whitaker, K. E., Bezanson, R., van Dokkum, P. G., et al. 
2017, ApJ, 838, 19, doi: 10.3847/1538-4357/aa6258 

Whitaker, K. E., Ashas, M., Illingworth, G., et al. 2019, 
ApJS, 244, 16, doi: 10.3847/1538-4365/ab3853 

Whitaker, K. E., Williams, C. C., Mowla, L., et al. 2021, 
Nature, 597, 485, doi: 10.1038/s41586-021-03806-7 

Wild, V., Charlot, S., Brinchmann, J., et al. 2011, MNRAS, 
417, 1760, doi: 10.1111/j.1365-2966.2011.19367.x 

Williams, C. C., Labbe, I., Spilker, J., et al. 2019, ApJ, 884, 
154, doi: 10.3847/1538-4357 /ab44aa 

Williams, C. C., Alberts, S., Ji, Z., et al. 2023a, arXiv 
e-prints, arXiv:2311.07483, 
doi: 10.48550/arXiv.2311.07483 

Williams, C. C., Tacchella, S., Maseda, M. V., et al. 2023b, 
arXiv e-prints, arXiv:2301.09780, 
doi: 10.48550/arXiv.2301.09780 


24 


Wuyts, S., Cox, T. J., Hayward, C. C., et al. 2010, ApJ, 
722, 1666, doi: 10.1088/0004-637X/722/2/1666 

Xiao, M. Y., Elbaz, D., Gómez-Guijarro, C., et al. 2023, 
A&A, 672, A18, doi: 10.1051/0004-6361/202245100 

Zavala, J. A., Casey, C. M., Manning, S. M., et al. 2021, 
ApJ, 909, 165, doi: 10.3847/1538-4357/abdb27 


Zhang, H., Primack, J. R., Faber, S. M., et al. 2019, 


MNRAS, 484, 5170, doi: 10.1093/mnras/stz339 


